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Abstract. We show that weak solutions of the relativistic Vlasov-Maxwell 
system preserve the total energy provided that the electromagnetic field is 
locally of bounded variation and, for any A > 0, the one-particle distribution 
function has a square integrable A-moment in the momentum variable. 



1. Introduction 

Consider an ensemble of relativistic charged particles that interact through their 
self-induced electromagnetic field. If collisions among the particles are so improba- 
ble that they can be neglected, then the ensemble can be modeled by the so-called 
relativistic Vlasov-Maxwell (RVM) system. At any given time t s]0, oo[, the RVM 
system is characterized by the one- particle distribution function / = f{t,x,p) with 
position X € K'^ and momentum p S R^. The self-induced electric and magnetic 
fields are denoted hy E = E{t, x) and B = B{t, x), respectively. Setting all physical 
constants to one, the model equations for a single particle species read 

df 

(1.1) -^+vVJ + iE + vxB)-S/pf^O 

(1.2) — - V X B = -Attj 
^ ' dt ■' 

(1.3) — + Vx^ = 

at 

(1.4) V-E^^Tip, V.B = 0, 

— 1/2 

where v := p (l -I- |p|^) denotes the relativistic velocity. The coupling of the 
Vlasov (|l.ip and Maxwell equations (|1.2p - p.4p is through the charge and current 
densities, which we denoted by p = p(t,x) and j — j{t,x) respectively. They are 
defined by 

(1.5) P fdp, j := vfdp. 

We define the Cauchy problem for the RVM system by (|l.ip - (|1.5p with initial data 

(1.6) /|t=o — fo, E^t^Q — Eo, -B|t=o = ^0, 

satisfying (|1.4p in the sense of distribution. It is not difficult to check that if (|1.4p 
holds at t = 0, then it will do so for all time in which the solution exist. Thus, the 
equations (jl.4p can be understood as a mere constraint on the initial data. 
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Now, define 

For T > 0, we say that (/, E, B) is a weak solution of the RVM system if 



i3 



(1.7) /GL°°([0,T[;4„nL-(R6)), E , B [L°" {%T[; L\R^))Y 

and the equations (|l.l[) - (|1.4p are satisfied in the sense of distributions. In particular, 
we say that the Vlasov equation p.ip is satisfied in the sense of distributions if for 
all ip e C°°([0,T] X R6) ^ith compact support in [0,T[xK6 

f{t, x,p) [dt^p + V ■ + K ■ Vpyj] (i, x,p)dtdxdp 



(1.8) =-/ fo{x,p)(p{0,x,p))dxdp 



We define analogous relations for the Maxwell equations (ll.2|l - (|1.4|) as well. The 
vector field K := E + v x B in ()1.8|) denotes the Lorentz force acting on a reference 
particle of velocity v. Notice that it satisfies Vp ■ K = 0. 

The global existence result for weak solutions in both relativistic and non- 
relativistic settings is due to DiPerna and Lions and can be found in [3]. In [7] 
this result is revisited. The uniqueness problem, on the other hand, remains un- 
solved. It is also unknown whether weak solutions preserve the total energy at least 
almost everywhere in time, cf. [31 Remark 4, p. 740]. It is known, however, that the 
energy is bounded at almost all t by its value at t — 0, namely 

eit) / Jl + \p\^f{t,x,p)dxdp 

■/R3xR3 

(1.9) +F- / \E{t,x)f + \B{t,x)\^dx < £{0). 

In the present note we show that if the electric and magnetic fields E and B are 
locally of bounded variation and, for any A > 0, the function 

(1.10) pxit,x) / \p\^ fit,x,p)dp 

is square integrable, then the relation (|1.9p is in fact an equality for almost all 
< t < T. Precisely, we prove the following result: 

Theorem 1. Let A > 0. Let fo G L^,„ n L^{R^) , Eo,Ba E [L^{M.^)]^ and 
denote by {f,E,B) a weak solution of the RVM system satisfying {f,E,B)\^^Q = 
{fo,Eo,Bo). IfE,B€ [Ll^{]Q,T[;BVioc{M.'^W and px as defined m [TW^ is m 
L,^^(]0,T[;L2(M3)), then the total energy defined by satisfies £{t) = £{Q) for 

almost all < t < T . 

The tools we use are basically those introduced by DiPerna and Lions in [3] 
to study renormalized solutions of transport equations. We shall also refer to [5], 
where applications to the Vlasov equation are given. We remark that the same 
result holds for the electromagnetic field in [-^^j^ocQO' '^[s ^/oc (^^))]^ since we have 
the (strict) inclusion W^'^{^) C BV{Q) for any open set fl C R^. For a detailed 
account on functions of bounded variation cf. [T] . We would Uke to include here the 
reference [S] , where the uniqueness of weak solutions for the Vlasov- Poisson system 
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has been obtained under the sole assumption that the spatial density is bounded. 
Similar results would be desirable for the more demanding Vlasov-Maxwell system. 

Formally, the law of the conservation of the total energy is derived as follows. 
Multiply the Maxwell equations (|1.2p and (|1.3|) by E and B respectively and inte- 
grate on M.^ to find that 

(1.11) / \E\' + \B\'dx^~ I j-Edx. 

57r dt ./b3 ./b3 



Multiply the Vlasov equation by y 1 + and integrate on M'' x to get 
(1.12) ^ f Jl + \p\^fdxdp ^ f j- Edx. 

Then the sum of (jl.lip and (|1.12p provide the desired result. 

As for weak solutions, we shall follow the same scheme. We find relations anal- 
ogous to (|l.lip and (|1.12p in sections [2] and [3] respectively. The difficulty is to 
overcome the lack of regularity and the need of justifying the operations taken for 
granted when the solutions are smooth. 

2. Energy balance for the Maxwell equation 

Here we show that if the current j is square integrable for almost all time, then 
the weak solution of the RVM system satisfies the energy balance associated to the 
Maxwell equations, i.e., the relation (|1.11[) . This result is reminiscent of the duality 
theorem for transport equations given by DiPerna and Lions in 4J. 

Lemma 1. Let (f,E,B) be a weak solution of the RVM system with initial data 
{fo,Eo,Bo). If j as defined in H^) is in [L°°(]0, T[; L'^{M.^))]^ , then 

(2.1) ^(||i?(t)|li. + WBml.J + 1^ ■ Edsdx = ^(lli^olli^ + ||i3o|li.) 
for almost all t G [0,T[. 

Proof. Let e > and let k e C^(R'^), k even, be a standard moUifier. Define the 
regularization kernel Ke ■— jtk(|-). Since mollification and distributional differenti- 
ation commute, i.e., (dxu) * = dx {u * k^), we can convolute (|1.2p and p.3p with 
Ke to obtain 

(2.2) _-l_VxB, = -47rje 

(2.3) -gf + VxE, = 0, 

where je ■= j * Ke and (i?e, Be) := {E, B) * k^. 

Consider the family of smooth cut-off functions (i>R :— 4'{^)i R > I where 
e C^(R^), > and (j) = 1 on the ball Bi C supp^ C -62- The smoothness of 
the fields Be and Ee with respect to x imply via (|2.2p and (|2.3p that dtEe, dtBe G 

[(]0,T[xR3)]3. ^^^^^ ^^^^^ ^ [wl^^^(]0,T[xR^) 
rule in Sobolev spaces, i.e., for almost all t g]0, T[ 
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and we can apply the chain 
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Therefore, we can multiply (12. 2|) and (12. 3|) by E^(f)R and i3e0K respectively, sum 
the resultant equations and integrate by parts to find that 



(2.4) = T- / / (^e X • V(/.fl - / f J,- E,<f>R. 

Let e — * 0. The terms on the left side converge as a consequence of the theorem 
of smooth approximations [51 Theorem 3, p. 196]. Also, the same theorem and the 
assumption made on the current j easily implies that J ■ E^ J j ■ E for almost 
all s e [0,T[. Thus, we may invoke the Lebesgue dominated convergence theorem 
and the convergence of the second term in the right side follows as well. Clearly, 
the same reasoning applies to the remaining term. Then, for almost all t G [0,T[ 



^ £^ [mf + \B{t)\') ^B-^j^^ {mf + is(o)i') H 

(2.5) = / / {ExB)-V^R- f f J- EcjiR. 

47r Jo JR3 Jo 



Finally, since for some constant Ct that does not depend on R 

I ||_B|| r^,2 , 



/ / {ExB)-V(j}R 
Jq Jm^ 



it is easy to check that (|2.ip follows from (|2.5p by letting R ^ oo. The proof of the 
lemma is complete. □ 



3. Energy Balance for the Vlasov Equation 
In this section we deduce the duality formula [4 resulting from the Vlasov equa- 
tion (jl.ip and (the identity) K -Vp^ 1 + \p\^ = v-E, which gives the energy balance 
associated to the Vlasov equation. Since we now face a nonlinear term in p.ip . we 
need to first prove the following lemma, a particular case of Lemma 3.5 in [5]. 

Lemma 2. Let k^-^ and k^^ be two regularization kernels defined on M.^ and 
respectively. Let (/, E, B) he a weak solution of the RVM system. Lf E, B £ 

[L^{]0,T[;BVioc(JS.^))]^, then there exist two sequences > 0, > 0, 0, 
£2 — > such that 

Va; ■ [v{k^^k^^ * /)] + Vp • [K{k^^k^^ * /)] 
(3.1) -(V, • [vf]) * K,^K,^ - (Vp • [Kf]) * 

converges to in L^{]Q,T[; Lj^^iR^ x R^)). 



Proof. First we omit the dependence in time and show the corresponding conver- 
gence on X R'^. Then we study the convergence on time as well. 
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Indeed, the compact support of the mohifiers and the divergence theorem allow 
us to rewrite (|3.ip as 

r(x,p) + /^(x,p) 

'J{x,p) - f{x -y,p- q)] [v{p) ~v{p-q)]- VyK^, {y)n^2 {Q)dydq 

[f{x,p) - f{x - y,p- q)][K{x,p) -K{x-y,p-q)] ■ V qK^^{q)n^^{y)dydq. 

In addition, since we have that 

K{x,p)- K{x-y,p-q) = E{x)~E{x-y)+v{p)x[B{x)~B{x-y)] 
(3.2) xB(a;-2/), 
we may decompose the second integral by 

where I^'^ involves the first two terms in the right side of p.2p and I^'^ involves 
the third term. Now, let i? > and define the set Br x Br =: fl C M.^ x such 
that Q + supp Kej^Kej C -B_R+i X Brj^i. In view of the assumptions of the lemma 

\\E{x) - E{x - y)\\Li(Ba) < W^-MMiB^^,) \yl \y\ < 

(similarly for B), where II Va;-B||_^^g^^^-j < oo denotes the norm of the measure 
WxE (resp. VxB), which coincides with the variation of E (resp. B) on the ball 
Bfl+i. Hence, since the relativistic velocity v € [C^(R'^)] satisfies |w| < 1, we find 
that for some positive constant Cr that depends on R 

lk'^'1lLi,^(0) ^ (ll^-^II.M(B„+i) + I|V-^IIa1(Sh+i)) 

(3-3) X /k2VgK,j) sup ||/(a;,p) - /(a; - q)||^^ . 



Similarly, we find the estimates 

(3-4) X /|e2VgKej) sup \\f{x,p) - fix -y,p- q)\\^2 (q) 

\J J lyl<^i,kl<^2 



and 



(3.5) X sup \!{x,p)~ j{x-y,p-q)\^2 ,m ■ 

Now, we have (J IeVkjI) < C, and we also have that 

sup \\f{x,p) - fix -y,p~ q)\\^2 {j« 0, as £1,62 ^0. 
Iy|<^i,l9l<s2 

Hence, we can choose two sequences £",£2 — * with £"/£2 = such that for 
some n sufficiently large the right-hand sides of p.4p and (|3.5p are less than 1/n. 
Therefore, since we also have / e L°°{^1), it follows that (|3.3p . (|3.4|) and (|3.5p go 
to zero as n — » 00, and so does p.ip in L;^^^(R'^ x M'^). 

Finally, we consider the dependence in time. The difficulty here seems to arise 
because the sequences e" and £2 may also depend on t. Otherwise we could just 
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invoke the Lebesgue dominated convergence theorem as p.3p and (|3.5p suggest. In 
particular, we must be careful with the estimate (13. 5|) . Nevertheless, if we keep 
track of the time dependence along the calculation, we find that 

II^^IIl;-^ ((o,T)xn) < C-r^ sup \\f{t,x,p)- f{t,x-y,p-q)\\^i,2 (^^o.T)xn) 

and we can reason as above. This concludes the proof of the lemma. □ 

We now turn to the energy balance (|1.12p associated to the Vlasov equation. 

Lemma 3. Let A > 0. In addition to the assumptions of Lemma\^ suppose that 
px as defined in ITJOj) is m L'^{]0,T[;L'^{R^)). Then 

(3.6) [ Jl + \p\^f{t)dxdp = / Jl + \p\^f{0)dxdp+ [ [ E ■ jdsdx 
Jr3xR3 Jr3xr3 Jo Jm^ 

for almost all t G [0,T[. 

Proof. (/, E, B) is a weak solution of the RVM system. Thus, as a straightforward 
consequence of Lemma [U there are two sequences e" > 0, > 0, e" — > 0, £2 ^ 
such that 

(3.7) 5tr + t;-V,/"+X-Vpr-r" 

converges to in {\Q,T[; Ll^^iM?' x M^)) as n ^ cx), where /" := k^^k^^ * / and 
r" is defined by (|3.ip . 

Consider a family of smooth cut-off hmctions =^ R ^ 1 where (f) G 



C^(R^), (/) > and = 1 on Si C supp0 C B2. If we multiply ([371) by ^ 1 + \pf(t)B. 
and integrate by parts, we find that 



1 + bi'/"(i)0fl - / \/i + \pff"mR 



E ■ vf'cj^R + / Jl + IpUnr'^ 
Jo Jk<5 

(3.8) + / / Jl + \pfrK ■ Vpct>R + f [ Jl + \pff"v ■ V.,(t,R. 

Jo JR6 JO JR6 

Here we have used the identity K ■ Vp^J 1 + = v ■ E. 

Let n —>■ 00. In doing so, we notice that the second term in the right-hand side 
vanishes as a consequence of Lemma[2] Also, the convergence of the two terms in the 
left-hand side and the last term in the right-hand side follow by a straightforward 
application of the theorem of smooth approximations. Thus, we are led to prove 
the convergence of the first and third terms in the right-hand side. 

Indeed, the reasoning done so far does not preclude us from writing (/(^ as the 
product of two suitable functions XR — x("r) ^^'^ Cr = C("r) where x G d^(M.^) 
and C e C^(M^). Hence, 

/ E-vif-fn^R < f \E\xrI CR\f-f'\ 

JR6 JR3 Jr3 

< c-fiiii^ii^, ii/-rii^. , 

a: X ,p 

which converges to zero as n — > 00. Then, a use of the Lebesgue dominated con- 
vergence theorem provide the convergence of the first term in the right-hand side. 
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Since we can do similarly with the remaining term, we find that as n — » oo, p.8|) 
converges to 

Ji + \p\\fmR = [ Ji + \p\^fmR + t f E- vf^R 

+ 11 \ll + \p?fK ■ Vp^R + f f Jl + \p\^fv ■ V^c^R. 

Jo JMB Jo JR8 

Finally, we let i? — > oo and show that the above equality converges to p.6p . The 
convergences of the term in the left and the first term in the right-hand side are 
straightforward, since for i = and for almost aU t > 0, f{t) e LJ.,;„(RS). Also, 
since 

ff \/i+\pffv-v.^R <§ r / \/i+\p\'f<^, 

Jo JR6 ^ Jo JR6 ^ 

the last term converges to zero as R —> oo. In order to obtain the convergence of 
the second term in the right, we first notice that for any A > 0, 



p{t,x) 



\p\<i 



f{t,x,p)dp 



f{t,x,p)dp 



\p\>i 



< oo. 



< y/4^\\fit,x)\\^,+px{t,x). 

Then, the hypothesis made in the lemma implies that 

(3.9) < VM3|I/WIIl2 +\\Pxit)h2 

As a result, and since \v\ < 1, we can easily verify that E ■ vf ^ L^Q^j r[xR^), so 
the Lebesgue theorem provides the expected convergence. Hence, we are only left 
to show that the third term in the right-hand side converges to zero. To this end, 
we first produce 



(3.10) 



/ / jl + \p\^fK-Vp(^R 

Jo Jmb 



< 



C 
R 



\E\ 



\B\) 



\p\f- 



R<\p\<2R 



To estimate the above inequality we observe that 



\p\f< 



R<\p\<2R 



< A < 1 

1 < A 



Thus, for any A > 0, there exists a constant C > independent of R such that the 
right-hand side of (|3.10p is less or equal than 



|pa(s)|Il2 ds < 



Therefore, l|3.10p converges to zero as i? — > cx) and the proof of the lemma is 
complete. □ 

4. Proof of Theorem [U 

Proof. Since |w| < 1, we have |j| < p. Then, in view of (|3.9p . we can combine 
Lemmas [1] and [3] to produce the equality for almost all < i < T claimed for 
Ol). □ 
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